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Direct-sequence  spread  spectrum  is  used  for  underwater  acoustic  communications  between  nodes,  at 
least  one  of  which  is  moving.  At-sea  data  show  that  the  phase  change  due  to  source  motion  is 
significant:  The  differential  phase  between  two  adjacent  symbols  is  often  larger  than  the  phase 
difference  between  symbols.  This  poses  a  challenge  to  phase-detection  based  receiver  algorithms 
when  the  source  or  receiver  is  moving.  A  pair  of  energy  detectors  that  are  insensitive  to  the  phase 
fluctuations  is  proposed,  whose  outputs  are  used  to  determine  the  relationship  between  adjacent 
symbols.  Good  performance  is  achieved  for  a  signal-to-noise  ratio  (SNR)  as  low  as  -10  dB  based 
on  at-sea  data.  While  the  method  can  be  applied  to  signaling  using  short  code  sequences,  the  focus 
in  this  paper  is  on  long  code  sequences  for  the  purpose  of  achieving  a  high  processing  gain  (at  the 
expense  of  a  low  data  rate),  so  that  communications  can  be  carried  out  at  a  low  input  SNR  to 
minimize  the  probability  of  detection  ( PD )  by  an  interceptor.  PD  is  calculated  for  a  typical  shallow 
water  environment  as  a  function  of  range  for  several  source  levels  assuming  a  broadband  energy 
detector  with  a  known  signal  bandwidth.  [DOl:  10.1121/1.2996329] 
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I.  INTRODUCTION 

Direct-sequence  spread-spectrum  (DSSS)  signaling  uses 
phase  coherent  signals  where  the  information  symbols  are 
multiplied  with  a  code  sequence,  commonly  known  as 
chips.1  The  signals  are  processed  at  the  receiver  using  the 
code  sequence  as  a  matched  filter  to  extract  the  information 
symbols.  '  Taking  advantage  of  the  processing  gain  derived 
from  the  matched  filter,  communications  can  be  carried  out 
at  low  signal  levels  (e.g.,  below  the  noise  level)  to  avoid 
detection  and  interception  by  an  unfriendly  party.  This  paper 
examines  the  probability  of  detection  and  counterdetection 
range  for  DSSS  signaling  in  an  underwater  acoustic  environ¬ 
ment  for  acoustic  communications  involving  mobile  plat¬ 
forms.  As  is  expected,  the  probability  of  detection,  as  well  as 
the  probability  of  interception,  is  very  much  dependent  on 
the  input  signal-to-noise  ratio  (SNR)  at  the  interceptor4'9  and 
the  recognition  differential  of  the  signal  waveform.  The  fact 
that  a  low  level  source  can  be  used  for  DSSS  communica¬ 
tions  and  that  the  signal  is  noiselike  presents  certain  advan¬ 
tages  from  the  point  of  view  of  low  probability  of  detection 
(LPD).  This  topic  will  be  analyzed  later  in  this  paper. 

The  problem  for  DSSS  communications  in  an  underwa¬ 
ter  acoustic  channel  is  the  multipath  arrivals,  which  create 
severe  interchip  and  intersymbol  interferences.  Decision 
feedback  equalizer  (DFE)  and  Rake  receiver  have  been 
adapted  for  DSSS  communications/'7  To  achieve  precise 
symbol  synchronization  and  channel  equalization,  high  SNR 
signals  are  required.8  Recently,  Yang  and  Yang91"  have  dem¬ 
onstrated  two  new  approaches  for  low  input-SNR  DSSS 
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communications  in  a  dynamic  ocean  environment  between 
fixed  nodes.  Minimal  (<1%)  bit  error  rate  (BER)  can  be 
achieved  for  in-band  SNR  as  low  as  -11  to  -14  dB.  Signal 
phase  distortion  due  to  signal  propagation  in  a  random  me¬ 
dium,  achievable  matched-filter  gain  (the  ratio  of  output  sig¬ 
nal  level  to  input  signal  level),  and  achievable  processing 
gain  (the  ratio  of  output  SNR  to  input  SNR)  were  measured 
from  data.  Performance  degradation  of  the  receiver  algorithm 
due  to  coarse  synchronization  at  low  input  SNRs,  signal  fad¬ 
ing,  and  imprecise  knowledge  of  the  channel  were  modeled 
and  compared  with  data. 

The  DSSS  signaling  method  uses  code  “orthogonality” 
to  minimize  interference  between  symbols.  The  code  or¬ 
thogonality  requires  that  the  code  sequence  is  almost  or¬ 
thogonal  to  any  of  the  cyclically  shifted  code  sequences. 
With  orthogonality,  the  matched-filter  output  yields  a  low 
sidelobe  level  and  thus  ensures  minimum  interference  in  a 
simple  environment  (dominated  by  a  single  path).  It  assures 
accurate  symbol  synchronization.  The  phase  of  the  matched- 
filter  output  can  then  be  used  to  determine  the  symbol  se¬ 
quence.  In  an  underwater  acoustic  channel  containing  many 
multipath  arrivals,  the  code  orthogonality  is  severely  de¬ 
graded  resulting  in  interchip  interference.  The  phase  is  often 
path  dependent  and  changing  rapidly  with  time.  At-sea  data 
showed  that  the  symbol  phase  error  estimated  from  the 
matched-filter  output  is  often  larger  than  the  phase  difference 
between  the  symbols,  resulting  in  unacceptable  bit  errors 
even  with  a  high  input  SNR.91"  The  challenge  is  to  develop 
a  receiver  algorithm,  patterned  after  the  simple  DSSS  signal 
processing  mentioned  above,  without  employing  a  DFE  to 
equalize  the  multipath  arrivals,  and  a  phase-locked  loop 
(PLL)  to  remove/compensate  the  phase  fluctuations  in  a  dy- 
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namic  ocean.  Toward  this  goal,  a  second  layer  of  correlator, 
which  cross-correlates  the  matched-filter  outputs  of  two  ad¬ 
jacent  symbols,  was  implemented  by  Yang  and  Yang,9’10  to 
estimate  their  differential  phase.  This  method  requires  only 
“coarse”  synchronization  at  the  symbol  level  and  is  appli¬ 
cable  to  low  input-SNR  communications  between  fixed 
nodes. 

In  practice,  underwater  acoustic  communications  often 
involve  mobile  platforms.  For  DSSS  signaling,  Doppler 
shift,  if  not  corrected  for,  destroys  the  orthogonality  of  the 
code  sequences  between  the  replica  signal  and  data,  and  be¬ 
tween  code  sequences  of  different  users.  In  addition,  the  sig¬ 
nal  coherence  time  is  significantly  reduced  due  to  the  source/ 
receiver  motion.11  The  reason  is  traced  to  the  rapidly 
fluctuating  (relative)  phases  of  the  multipaths.  As  a  symbol 
phase  is  determined  by  interference  between  the  multipaths, 
medium  fluctuation  causes  a  symbol  phase  to  change  with 
time.  For  a  moving  source,  the  individual  paths  encounter,  in 
addition,  a  different  phase  change  (ray  path-length  change) 
due  to  source  changing  range.  A  higher  rate  of  phase  fluctua¬ 
tion  is  thus  expected,  resulting  in  a  shorter  coherence  time 
for  a  moving  source  than  for  a  fixed  source.  One  finds9'10  for 
a  fixed  source/receiver  pair  that  while  the  phase  of  individual 
symbols  changes  rapidly  with  time,  the  differential  phase 
between  the  symbols  is  relatively  small.  In  contrast,  the 
analysis  below  will  show  that,  for  the  moving-source  data, 
after  proper  Doppler  compensation,  the  differential  phase  be¬ 
tween  symbols  changes  rapidly  with  time,  and  is  often  larger 
than  the  phase  difference  between  true  symbols  even  with  a 
high  input  SNR.  Hence,  while  differential  phase  shift  keying 
(DPSK)  can  be  used  to  decode  the  signal  for  the  fixed-source 
data,  this  method  (based  on  the  differential  phase)  is  ex¬ 
pected  to  encounter  high  bit  errors  when  applied  to  the 
moving-source  data. 

For  DSSS  signaling  involving  a  mobile  platform,  a  new 
method,  which  is  insensitive  to  the  phase  fluctuation,  is  de¬ 
veloped  in  this  paper  and  applied  to  at-sea  data.  This  method 
can  also  be  applied  to  fixed-source  data.  It  will  be  shown 
below  that  the  method,  applied  to  the  fixed-source  data  re¬ 
ported  previously,  yields  a  slightly  degraded,  but  still  satis¬ 
factory,  performance  compared  with  the  method  mentioned 
above. 

Like  the  cross-correlation  method,  the  new  method  can 
be  applied  to  DSSS  signals  with  short  or  long  code  se¬ 
quences  so  long  as  the  code  length  is  longer  than  the  multi- 
path  spread.  Similar  to  the  previous  work,9’10  the  focus  of 
this  paper  is  on  communications  with  a  long  code  sequence 
for  the  purpose  of  LPD.  A  long  sequence  yields  a  high  pro¬ 
cessing  gain,  which  allows  communications  using  a  low 
source  level  (SL),  thus  minimizing  the  probability  of  detec¬ 
tion  by  a  hostile  party.9,10  For  a  given  bandwidth,  longer 
sequences  imply  lower  data  rates.  One  can  increase  the  data 
rate  through  other  means,  such  as  spreading  the  data  over 
multiple  orthogonal  codes  at  the  price  of  a  slightly  higher 
SL.  That  topic  is  beyond  the  scope  of  this  paper. 

This  paper  is  organized  as  follows:  In  Sec.  II  the  char¬ 
acteristics  of  a  signal  emitted  from  a  moving  source  are  ana¬ 
lyzed,  including  the  fluctuation  of  the  Doppler  shift  found  in 
the  data.  In  Sec.  Ill,  the  phase  insensitive  method  is  pre- 


km 

FIG.  1 .  (Color  online)  Source-receiver  geometry  including  the  towed  source 
run. 

sented,  tested  against  the  fixed-source  data,  and  then  applied 
to  the  moving-source  data.  In  Sec.  IV,  the  probability  of  de¬ 
tection  and  the  probability  of  false  alarm  are  analyzed  for  an 
energy  detector  assuming  the  interceptor  has  no  knowledge 
of  the  signal  waveform  and  signal  parameter.  The  ability  to 
communicate  with  a  lower  SL  using  the  DSSS  signaling  is 
shown  to  significantly  reduce  the  probability  of  detection  by 
an  interceptor. 

Regarding  LPD  underwater  acoustic  communications,  a 
simple  system  concept  may  be  helpful  to  motivate  and  guide 
the  analysis  below.  A  potential  scenario  for  LPD  communi¬ 
cations  will  be  cooperative  operation  between  autonomous 
underwater  vehicles  (AUVs),  which  stay  together  within  a 
certain  range  from  each  other  and  use  LPD  acoustic  commu¬ 
nications  to  avoid  detection  by  a  nearby  interceptor.  From 
the  system  point  of  view,  the  SL  needs  to  be  high  enough  to 
assure  reliable  communications  with  neighboring  AUVs,  and 
yet  low  enough  so  not  to  alert  the  interceptor  of  their  pres¬ 
ence.  The  question  is  at  what  range  will  the  detection  (by  the 
interceptor)  become  unavoidable,  i.e.,  the  counterdetection 
range.  Counterdetection  range  is  calculated  in  Sec.  IV  for  a 
typical  shallow  water  environment  assuming  certain  system 
parameters. 

II.  SIGNAL  CHARACTERISTICS 

The  data  were  taken  during  the  TREX04  (time  reversal 
experiment)  conducted  by  the  Naval  Research  Laboratory  in 
April  2004,  off  the  coast  of  New  Jersey,  southwest  of  the 
Hudson  Canyon.  The  sound  speed  profile  presented  a  mini¬ 
mum  in  the  middle  of  the  water  column  (see  Fig.  1,  Ref.  10). 
Acoustic  communication  data  were  transmitted  from  a  fixed 
source  to  a  fixed-receiver  array  at  the  range  of  3.4  km,  and 
from  a  towed  source  to  the  fixed  receiver  as  shown  in  Fig.  1. 
The  water  depth  in  the  experimental  area  is  about  70  m.  The 
fixed  source  and  receivers  were  located  at  about  35  m  depth. 
The  towed  source  was  at  a  depth  of  approximately  20  m.  The 
vertical  array  has  an  aperture  of  approximately  2  m  and  con¬ 
tains  eight  hydrophones  with  nonuniform  spacing.  Only  one 
receiver  is  used  for  communications  from  the  fixed  and 
towed  sources. 
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FIG.  2.  (Color  online)  Estimation  of  the  mean  Doppler  shift  over  a  packet  as 
shift  within  a  packet  based  on  ///-sequence  ambiguity  function  (b). 

The  DSSS  signals  were  centered  at  17  kHz  and  had  a 
bandwidth  of  4  kHz.  The  transmitted  symbols  were  spread 
with  (multiplied  by)  an  m-sequence  containing  511  chips 
yielding  a  theoretical  processing  gain  of  27  dB  at  the  re¬ 
ceiver.  The  sequence  of  chips  coded  with  the  message  data 
was  then  transmitted  using  binary  phase-shift-keying  modu¬ 
lation.  Each  symbol  has  a  duration  of  127.8  ms,  yielding  a 
data  rate  of  ~8  bits/s.  So  long  as  the  symbol  duration  is 
longer  than  the  multipath  spread,  the  interference  between 
adjacent  symbols  is  minimal.  [For  example,  see  Fig.  1(b)  of 
Ref.  10.]  Each  packet  contains  200  symbols  lasting  approxi¬ 
mately  25.55  s.  A  total  of  111  packets  were  transmitted  over 
approximately  an  hour.  The  signal  properties  are  discussed  in 
this  section.  BERs  are  reported  in  the  next  section. 

In  this  experiment,  some  auxiliary  signals,  such  as  pilot 
tones  and  linear  frequency  modulation  (LFM)  signals,  were 
added  to  each  packet  for  independent  probing  of  the  channel 
properties.  Pilot  tones  signals  shall  not  be  used  for  DSSS 
communications  at  a  low  SNR.  Data  acquisition  and  symbol 
synchronization  can  be  done  by  matched  filtering  the  data 
with  the  LFM  or  m-sequence  replicas.9'10 

A.  Doppler  shift 

A  key  feature  of  the  moving-source  data  is  the  Doppler 
shift/spread.  The  Doppler  shift  is  path  dependent  as  multi¬ 
paths  arrive  with  different  depression/elevation  (D/E)  angles. 
At  long  ranges,  the  multipaths  are  confined  to  shallow 
(<10°)  D/E  angles  as  deep  D/E  angle  paths  are  attenuated 
by  the  bottom.  Under  these  conditions,  one  can  assume  the 
same  Doppler  shift  for  all  the  paths.  At  short  ranges,  both 
shallow  and  deep  rays  exist.  As  a  result,  the  path-dependent 
Doppler  shifts  create  a  (large)  Doppler  spread.  Note  that  an 
individual  path  may  also  have  a  Doppler  spread  due  to  scat¬ 
tering  with  ocean  medium  and/or  surface.  The  path- 
dependent  Doppler  spread,  most  notable  at  short  ranges,  can 
be  (much)  larger  than  the  medium/boundary  induced  Dop¬ 
pler  spread. 
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function  of  packet  number  using  different  signals  (a).  Instantaneous  Doppler 

The  first  task  in  the  data  analysis  is  Doppler  estimation 
as  the  data  have  to  be  corrected  for  the  Doppler  shift.  The 
temporal  variation  of  the  Doppler  shift  on  the  scale  of  the 
order  of  a  minute  can  be  obtained  by  measuring  the  average 
Doppler  shift  for  each  packet.  Three  methods  are  used  to 
estimate  the  average  Doppler  shift.  The  first  method  deter¬ 
mines  the  average  Doppler  shift  from  the  pilot  tones,  which 
were  placed  above  and  below  the  DSSS  signal  band.  The 
result  is  shown  as  the  dotted  line  in  Fig.  2(a).  The  second 
method  determines  the  average  Doppler  shift  from  the  time 
separation  between  the  two  LFM  signals  placed  before  and 
after  the  data  packet  measured  against  that  of  the  transmitted 
signal.  The  measured  Doppler  shift  is  shown  as  the  dashed 
line  in  Fig.  2(a).  The  third  method  measures  the  Doppler 
shift  from  each  m-sequence.  The  input  data  are  bandpass 
filtered  and  shifted  to  the  base  band.  The  base  band  data  are 
synchronized  and  partitioned  into  (overlaying)  blocks,  each 
containing  one  m-sequence.  A  wideband  ambiguity  function 
is  created  for  each  block  of  data  by  correlating  the  data  with 
a  Doppler-shifted  transmitted  m-sequence  signal.  The  output 
is  an  impulse  response  as  a  function  of  Doppler  shift  and 
delay  time,  referred  to  as  the  ambiguity  function.  The  peak  of 
the  ambiguity  surface  determines  the  Doppler  shift  for  the 
data  block.  The  measured  Doppler  shift  as  a  function  of  the 
block  (or  sequence)  number  is  shown  in  Fig.  2(b)  using  one 
packet  of  data  as  an  example.  The  Doppler  shift  is  averaged 
for  each  packet,  and  the  result  is  shown  in  Fig.  2(a)  by  the 
solid  line  as  a  function  of  packet  number. 

Based  on  the  source  track  in  Fig.  1,  the  source  moved 
out  initially  by  an  opening  range.  It  turned  around  at  a  range 
about  5  km  and  moved  toward  the  receiver;  it  reached  a  clos¬ 
est  point  of  approach  near  the  end  of  the  run.  [The  source- 
receiver  range  is  shown  in  Fig.  3(a)  as  a  function  of  the 
packet  number.]  One  thus  expects  initially  a  negative  Dop¬ 
pler  shift,  which  should  change  into  a  positive  Doppler  shift 
and  then  crosses  over  into  a  negative  Doppler  shift  again. 
This  is  in  accord  with  the  measured  Doppler  shift.  One  finds 
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Symbol  No 

FIG.  4.  (Color  online)  Signal  amplitude  as  a  function  of  time  within  a 
packet  for  three  selected  receivers. 


FIG.  3.  (Color  online)  Source-receiver  range  as  a  function  of  time  (packet 
number)  (a).  Signal  coherence  time  as  a  function  of  time  (packet  number) 
(b). 

that  the  Doppler  shift  measured  from  the  LFM  and 
m-sequence  signals  agree  relatively  well  with  each  other.  The 
measurement  using  the  pilot  tones  is  somewhat  (a  few  hertz) 
lower  than  using  LFM  and  m-sequence.  The  reason  is  the 
ambiguity  in  the  spectral  determination  of  the  Doppler  shift. 
Using  the  spectral  analysis,  an  integration  time  of  2  s  is  re¬ 
quired,  for  example,  to  achieve  a  Doppler  resolution  of 
0.5  Hz.  However,  during  the  2  s  period,  the  Doppler  shift 
changes  with  time  as  shown  next. 

Figure  2(b)  shows  the  Doppler  shift  within  a  packet. 
One  may  expect  the  Doppler  shift  to  vary  slowly  with  time 
[based  on  Fig.  2(a)],  and  be  surprised  that  the  Doppler  shift 
varies  by  3-5  Hz  with  a  period  of  4-7  s.  This  variation  is 
attributed  to  the  heave  motion  of  the  ship,  which  was  ob¬ 
served  during  the  experiment,  and  expected  to  cause  the  tow 
body  to  speed  up  and  down  with  the  ship  motion.  For  data 
processing,  what  is  significant  is  that  the  Doppler  shift  varies 
from  symbol  to  symbol  [Fig.  2(b)];  each  symbol  is  —128  ms 
long.  This  requires  the  Doppler  shift  to  be  estimated  for  ev¬ 
ery  symbol,  which  increases  significantly  the  data  processing 
load. 

6.  Temporal  coherence:  Signal  coherence  time 

The  Doppler  corrected  impulse  responses  are  used  to 
estimate  the  temporal  coherence  of  the  m-sequence  signals  as 
defined  by 

,  ^  [p*(t)p(t+  t)] 

P  +  T)p(t  +  r)] 

where  p(t)  and  p{t+  t)  are  the  impulse  responses  of  the  two 
consecutive  m-sequences  separated  by  a  lag  time  r,  t  is  the 
sample  time  of  the  impulse  response.  The  square  bracket 
[af>]  denotes  the  maximum  absolute  value  of  the  correlation 
between  a(t )  and  b(t).  (For  the  measurement  of  the  impulse 
response  and  temporal  coherence  using  m-sequences,  the 
reader  is  referred  to  Ref.  13.)  The  signal  coherence  time  r0  8 


is  defined  as  the  lag  time  when  the  coherence  value  drops 
from  1  to  0.8.  Temporal  coherence  of  the  m-sequence  signal 
is  measured  for  each  packet  and  the  measured  signal  coher¬ 
ence  time  is  shown  as  a  function  of  packet  number  in  Fig. 
3(b). 

One  finds  in  Fig.  3(b)  that  the  signal  coherence  time  at 
close  ranges  (<1  km)  is  —0.1  s  compared  with  the  signal 
coherence  time  at  3  km,  which  is  of  the  order  of  0.5  s.  Note 
that  the  signal  coherence  time  is  inversely  proportional  to  the 
signal  Doppler  spread.  At  close  ranges,  the  multipaths  have 
different  Doppler  shifts,  resulting  in  a  decreased  coherence 
time.  Later  in  Sec.  Ill,  it  will  be  shown  that  the  “large” 
Doppler  spread  presents  a  problem  for  Doppler  compensa¬ 
tion  and  is  responsible  for  degraded  BER  performance  at 
close  ranges. 


C.  Signal-fading  statistics 


Figure  4  shows  an  example  of  the  time  variation  of  the 
symbol  amplitude  using  one  packet  of  data.  The  amplitude 
variation  is  shown  for  three  receivers  spaced  by  approxi¬ 
mately  0.3  m,  as  well  as  the  amplitude  of  the  linear  sum  of 
the  pressure  fields  on  all  eight  receivers.  Fading  by  as  much 
as  10-20  dB  is  seen  in  the  data.  Since  detection  performance 
is  affected  by  signal  fading,  one  needs  to  examine  the  signal¬ 
fading  statistics  for  a  moving  source  (on  a  single  receiver). 
For  amplitude  statistics,  the  amplitude  is  estimated  from  the 
square  root  of  the  symbol  energy,  conventionally  denoted  as 
Eh.  Figure  5(a)  shows  the  symbol  amplitude  distribution  at  a 
fixed  range  of  3.4  km  (Fig.  10  of  Ref.  10).  It  was  found 
consistent  with  a  log-normal  distribution 

(Ins- 


p(s): 


1 


V2iro 


exp 


2(77 


(2) 


where  s  is  the  signal  amplitude;  s>0.  The  mean  and  vari¬ 
ance  of  the  signal  amplitude  are  expressed  in  terms  of  two 
parameters,  fix  and  rrt,  as  follows: 

ms  =  exp(yu,s  +  ct2/2)  , 
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per  symbol 


FIG.  5.  (Color  online)  Probability  distribution  of  signal  amplitude  for  fixed-source  data  compared  with  the  log-normal  distribution  (a)  and  histogram 
distribution  of  signal  amplitude  from  the  moving-source  data  (b). 


v;  =  exp(2 fis  +  cr3)[exp(cr3)  -  1]. 

The  data  were  fitted  in  Fig.  5(a)  using  fis=- 1.6449  and  a2 
=  0.0304;  the  corresponding  mean  signal  level  and  its  vari¬ 
ance  are  «is= 0.196  and  u3=1.186X  10-3. 

One  can  generalize  Eq.  (2)  to  model  the  data  from  a 
moving  source,  such  as  communication  data  from  an  autono¬ 
mous  underwater  vehicle.  As  the  range  changes  from  rx  to  r2, 
the  amplitude-fading  statistics  will  be  a  weighted  average  as 

rm  2 

p(s)=  p(ms)ps(s\ms)dms,  (3) 

J  W] 

where  mx  and  m2  are  the  mean  signal  level  at  ranges  rx  and 
r2,  respectively,  p{ms )  is  the  probability  function  for  the 
mean  signal  level  over  the  range  track,  and  ps(s\ms)  is  the 
amplitude-fading  distribution  given  in  Eq.  (2).  The  resulting 
amplitude-fading  statistics  is  not  expected  to  be  log-normal 
unless  m  is  nearly  constant. 

Figure  5(b)  shows  the  histogram  distribution  of  the  sym¬ 
bol  amplitude  determined  from  the  packets  transmitted  when 
the  source  moves  from  5  to  0.5  km.  Compared  with  Fig. 
5(a),  one  finds  that  the  two  distributions  are  similar.  One 
observes  that  the  mean  signal  energy  of  the  packets  fluctu¬ 
ates  by  ±3  dB  over  the  track.  This  small  variation  (limited 
range  change)  is  the  reason  for  the  qualitatively  similar  dis¬ 
tributions  between  Figs.  5(a)  and  5(b). 

III.  PHASE  ERROR,  TRANSITION  DETECTOR,  AND 
BER  ANALYSIS 

A.  Symbol  phase  error 

Figure  6(a)  shows  the  differential-phase  error  between 
symbols  (the  error  of  the  differential  phase  of  symbols  be¬ 
tween  data  and  transmitted  signal)  determined  from  the  fixed 
source/receiver  data  [repeated  from  Fig.  6(b)  of  Ref.  10].  As 
discussed  above,  the  differential  phase  is  determined  from 
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the  peak  of  the  cross  correlation  of  the  matched-filter  outputs 
of  two  adjacent  symbols.  The  cross  correlation  of  the 
matched-filter  outputs  yields  (h.  hi+i)S .  Si+X  plus  terms  con- 

•  •  •  10  i  i 

taining  noise,  which  will  be  ignored  for  the  moment,  where 
hj  is  the  channel  impulse  response  for  the  ith  symbol,  de¬ 
noted  by  S; .  Assuming  that  the  channel  has  not  changed 
much  within  the  time  period  of  two  symbols,  i.e.,  h.hi+l~I, 
one  can  estimate  the  differential  phase  between  the  two  sym¬ 
bols  from  the  cross-correlation  result.  This  is  a  simple  way  to 
estimate  the  symbol  differential  phase  without  requiring  a 
channel  equalizer.  The  differential-phase  error  [shown  in  Fig. 
6(a)]  is  obtained  by  subtracting  the  true  differential  phase 
(between  the  two  symbols)  from  the  measured  differential 
phase.  The  small  differential-phase  error  (smaller  than  the 
phase  difference  between  symbols,  i.e.,  ±90°)  suggests  that 
for  the  fixed-source  data,  the  symbol  phase  error  can  be  cor¬ 
rected  using  a  PLL.  For  practical  applications  where  process¬ 
ing  power  is  limited,  it  is  much  simpler  to  use  DPSK  signal¬ 
ing  without  requiring  the  PLL.  For  the  data  shown  in  Fig. 
6(a),  there  is  no  bit  error  since  the  phase  error  is  <  ±90°. 

When  either  of  the  source  and  receiver  is  moving,  one 
needs  to  estimate  the  Doppler  shift  for  each  symbol  (as  dis¬ 
cussed  above)  and  apply  Doppler  correction  (including  both 
carrier  phase  offset  and  symbol  dilatation)  to  the  data.  A 
block  diagram  for  the  processing  is  shown  in  Fig.  7.  For 
details  see  Sec.  Ill  C.  Applying  this  method  to  the  moving- 
source  data,  one  can  deduce  the  differential-phase  error  of 
symbols  as  shown  in  Fig.  6(b).  It  shows  that  a  high  percent¬ 
age  of  the  symbols  are  in  error  for  the  moving-source  data. 
As  explained  above,  the  symbol  phase  is  rapidly  changing 
for  a  moving  source  since  the  phase  of  each  individual  path 
is  varying  with  time.  The  interference  among  the  time- 
varying  multipath  arrivals  results  in  a  fast  changing  differen¬ 
tial  phase  between  symbols.  To  mitigate  the  fast  changing 
symbol  phase,  an  energy  detector  that  is  insensitive  to  the 
symbol  phase  fluctuations  is  proposed  below. 
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FIG.  6.  (Color  online)  Differential  symbol  phase  error  as  function  of  symbol  number:  fixed-source  data  (a)  and  moving-source  data  (b).  Fixed-source  data  are 
reproduced  from  Ref.  10.  Moving-source  data  points  (X)  are  connected  by  lines  to  show  the  magnitude  of  the  phase  variation. 


B.  Transition  detectors 

In  a  time-invariant  environment  without  multipaths,  one 
could  determine  the  next  symbol  from  the  current  symbol 
using  a  pair  of  code  sequences  proposed  previously,14 

T+  =  [C,C]  and  T_=[C,-C],  (4) 

where  C  is  the  code  sequence  expressed  as  a  row  vector.  The 
data,  after  proper  synchronization  and  Doppler  correction, 
are  divided  into  blocks.  The  data  blocks  containing  two  sym¬ 
bols  can  be  expressed  as 

d  =  [y„Ce'e",5„+1Ce,6l"+l]  +  N,  (5) 

where  s„  is  the  nth  symbol,  sn=  ±  1,  0n  is  the  random  symbol 
phase,  and  N  denotes  the  noise;  in  Eq.  (5)  the  signal  ampli¬ 
tude  is  normalized  to  1.  The  matched-filter  output  between 
Eqs.  (4)  and  (5)  yields 

T+  ■  d  =  ( s„e'e»  +  sn+le‘e^)p  +  pN, 
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FIG.  7.  Block  diagram  of  the  cross-correlation  method  for  DPSK  signals. 


T_-d=  ( sne'e « -  s„+le'e"+l)p+  pN,  (6) 

where  p  is  the  autocorrelation  function  of  the  code  sequence, 
p=OC,  and  pN  denotes  the  correlation  of  the  code  sequence 
with  the  noise,  pN=C'N.  Since  p  >  pN  assuming  a  high  pro¬ 
cessing  gain,  the  noise  contribution  is  negligible  except  at 
very  low  SNRs.  For  simplicity  we  shall  drop  the  noise  in  the 
discussions  throughout  this  paper. 

In  a  channel  without  or  with  minimal  differential-phase 
error,  i.e.,  6n=9n+\ ,  one  finds  max|7’+,c/|2>max|7’_,£/|2,  if 
s„s„+l  =  1  (adjacent  symbols  remain  the  same,  no  transition), 
and  max|7’_-d|2>max|7’+'<af|2,  if  s„s„+ 1  =  - 1  (adjacent  symbol 
not  the  same).  Thus  by  comparing  the  outputs  of  the  matched 
filters,  Eq.  (6),  one  can  determine  the  next  symbol  based  on 
the  current  symbol.  However,  one  finds  that  this  method  may 
not  work  in  a  channel  that  contains  a  significant  differential- 
phase  error  since  the  output  squared  of  the  matched  filter 
[Eq.  (6)]  is  a  function  of  the  symbol  differential  phase, 
which,  when  in  error,  can  cause  incorrect  symbol  decisions. 

To  process  the  moving-source  data,  which  contain  a  sig¬ 
nificant  differential-phase  error,  a  new  pair  of  transition  de¬ 
tector  is  proposed,  which  is  insensitive  to  the  differential- 
phase  error  and  therefore  robust  for  symbol  estimation/ 
decision.  Let  Cj  be  the  first  half  of  the  code  sequence,  and  C2 
be  the  second  half  of  the  code  sequence,  i.e.,  C=[C[,C2]. 
The  new  transition  detectors  are  given  by 

Cp=[C2,Cl ]  and  CN=  [C2,~  C,].  (7) 

To  see  how  the  new  transition  detector  works,  one  compares 
with  the  (conventional)  matched-filter  operation.  In  a  multi- 
path  time-varying  environment,  the  data  of  Eq.  (5)  can  be 
expressed  in  four  half-blocks  as 

d=  [snhn  Chsnhn+\i2C2,sn+fhn+[C  hsn+fhn+ut  /2C2] ,  (8) 

where  hn  is  the  impulse  response  at  f„,  and  /z„+1/2  is  the  im¬ 
pulse  response  at  tn+rl 2;  t  is  the  symbol  duration.  The 
matched  filter  (using  the  code  sequence  C)  yields 
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FIG.  8.  (Color  online)  BER  (a)  and  average  BER  (b)  as  a  function  of  input  SNR  using  the  cross-correlation  method  for  fixed-source  data.  The  solid  line  is 
a  model  prediction.  See  Ref.  10.  Zero  BER  is  represented  by  10-6  on  the  logarithmic  scale. 


C  ■  d~  \snhnp,sn+ihn+ip\,  (9) 

where  one  approximates  hn,  hn+l/ 2  by  h„;  the  conventional 
matched  filter  assumes  that  the  channel  is  time  invariant  dur¬ 
ing  a  symbol  period.  The  correlation  of  the  matched-filter 
outputs  between  the  two  symbols  yields 

(s*nsn+l)(h*n  ■  hn+1)(p .  p )  =  (s*sn+i)\h*hn+i\el<l>n+l{p  ■  p), 

(10) 

where  <b  is  the  differential  phase  between  the  symbol  due  to 
the  random  medium,  which  causes  the  differential-phase  er¬ 
ror. 

The  new  transition  filter  works  by  the  same  principle  as 
the  (conventional)  matched  filter,  Eq.  (9).  One  finds 

C±-d=hn+mC±  •  [.s„C2,VnC|],  (11) 


where  d  contains  two  half-blocks  of  data,  i.e.,  d 
=  \sJ1n+mC2^sn+i^ln+iCi],  again  assuming  that  the  channel  is 
invariant  during  one  symbol  block  d.  One  can  determine 
whether  the  next  symbol  is  the  same  as  (or  opposite  to)  the 
current  symbol  by  comparing  the  energy  of  the  above  tran¬ 
sition  detector.  For  example,  if  sn+l=sn,  one  has 


max|C±  •  d\2 


sA+i/2|2(max|C±  •  C+ 1)2 


(12) 


where  the  autocorrelation  function  of  C+  is  the  same  as  the 
autocorrelation  function  of  C  itself;  it  has  a  peak  value  equal 
to  the  length  of  the  code  sequence,  M,  and  a  sidelobe  level 
equal  to  1,  and  the  correlation  of  C_  against  C+  has  a  peak 
value  (i  much  less  than  M;  the  actual  value  depends  on  the 
sequence  used.  Thus  the  C+  matched  filter  will  yield  a  higher 
peak  energy  than  the  C_  matched  filter.  Conversely,  if  the 
peak  energy  of  the  C+  matched-filter  output  is  higher  than 


that  of  the  C_  matched  filter,  it  would  suggest  that  sn+\  =sn. 
Likewise,  if  sn+\=-sn,  one  has 


max|C±  ■  d\2 


-  Mh+I/zI 

=  |/bi+l/2|2 


2(max|C±  •  C_|)2 

V 

M 2  ’ 


(13) 


where  the  correlation  of  C_  against  itself  has  a  peak  value 
equal  to  the  length  of  the  code  sequence  and  a  sidelobe  level 
that  varies  depending  on  the  sequence;  normally  close  to  1. 
The  correlation  of  C+  against  C_  has  a  peak  value  f3  much 
less  than  M  as  mentioned  above.  Since  M>(3,  the  C_ 
matched  filter  will  yield  a  higher  peak  energy  than  the  C+ 
matched  filter.  Conversely,  if  the  peak  energy  of  the  C_ 
matched-filter  output  is  higher  than  that  of  the  C+  matched 
filter,  it  would  suggest  that  sn+\=-sn. 

Note  that  the  outputs  of  Eqs.  (12)  and  (13)  are  real  num¬ 
bers.  The  proposed  transition  C±  detector  is  thus  not  affected 
by  the  problem  of  rapid  phase  variations  and  can  be  applied 
to  the  moving-source  data.  The  results  will  be  shown  in  Secs. 
Ill  C  and  III  D.  The  transition  detector  can  also  be  applied  to 
the  fixed-source  data.  We  will  first  analyze  the  fixed-source 
data  using  the  transition  detector  to  compare  its  performance 
and  with  the  cross-correlation  method  of  Ref.  10. 

The  fixed-source  data  contained  1160  packets  with  an 
input  SNR  varying  from  -15  to  23  dB.  The  BER  result  by 
cross-correlating  the  matched-filter  outputs  of  two  adjacent 
symbols  is  shown  in  Fig.  8  (Fig.  8  of  Ref.  10)  assuming 
DPSK  signaling.  Figure  8(a)  shows  the  BER  as  a  function  of 
the  packet  input  SNR  (the  input  SNR  averaged  over  each 
packet)  and  Fig.  8(b)  shows  the  BER  averaged  over  packets 
with  approximately  the  same  input  SNR,  plotted  as  a  func¬ 
tion  of  the  input  SNR.  Also  shown  in  Fig.  8(a)  is  the  mod¬ 
eled  BER  for  the  DPSK  signal.  Using  the  transition  filter,  the 
BER  result  is  shown  in  Fig.  9(a)  as  a  function  of  the  packet 
input  SNR,  and  the  average  BER  is  shown  in  Fig.  9(b)  as  a 
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FIG.  9.  (Color  online)  BER  (a)  and  average  BER  (b)  as  a  function  of  input  SNR  using  the  transition  detector  method  applied  to  the  fixed-source  data.  X 
denotes  the  data  of  Fig.  8(b)  using  the  cross-correlation  method.  The  solid  line  is  a  model  prediction.  Zero  BER  is  represented  by  10-6  on  the  logarithmic  scale. 


function  of  the  input  SNR;  both  denoted  by  Also  shown 
in  Fig.  9(b)  is  the  average  BER  repeated  from  Fig.  8(b), 
shown  by  “x.”  One  finds  that  the  BER  results  using  the  tran¬ 
sition  detector  are  slightly  higher  than  that  using  the  cross 
correlation  of  the  matched-filter  outputs.  For  the  latter 
method,  as  pointed  out  in  Ref.  10,  the  phase  estimation  can 
be  improved  under  certain  channel  conditions  by  synchroniz¬ 
ing  the  cross-correlation  output  based  on  symbols  with 
higher  SNRs.  The  BER  in  the  actual  data  can  be  better  than 
expected  (modeled)  whenever  the  symbols  happen  to  be 
properly  synchronized.  This  feature  is  not  available  to  the 
transition  detector  (an  energy  detector)  for  LPD  applications. 

C.  Moving-source  data  results:  High  SNR 

The  moving-source  data  require  initial  data  processing 
involving  Doppler  estimation  and  correction.  The  block  dia¬ 
gram  of  data  processing  is  shown  in  Fig.  10.  Doppler  esti¬ 
mation  is  done  in  the  base  band  with  two  samples  with  frac¬ 
tional  sampling  (two  samples  per  chip).  After  coarse 

Coarse  Synchronization  Received  data  partitioned  into  blocks 


FIG.  10.  (Color  online)  Block  diagram  of  the  transition  detector  method.  Cl 
and  C2  are  the  first  and  second  halves  of  the  code  sequence  as  described  in 
Sec.  IH  B. 


synchronization,  the  data  are  divided  into  overlapping  blocks 
with  a  length  larger  than  the  code  length  to  include  the  syn¬ 
chronization  error.  Doppler  shift  is  estimated  for  each  block 
(symbol)  using  the  wideband  ambiguity  functions  as  de¬ 
scribed  above.  The  signal  phase  shift  (due  to  Doppler  shift) 
is  corrected  for  each  block  of  data  based  on  the  estimated 
Doppler  shift.  The  data  are  resampled  based  on  the  estimated 
Doppler  shift  to  correct  for  time  dilation  or  compression  and 
resynchronized.  The  symbols  are  detected  by  applying  the 
transition  filter  C±  to  the  data  to  determine  if  the  two  adja¬ 
cent  symbols  are  the  same  or  not.  In  actual  implementations, 
one  finds  other  ways  to  simplify  the  processing.  For  ex¬ 
ample,  it  is  simpler  to  Doppler-shift  the  transition  filter  C± 
and  apply  the  filter  to  the  data  according  to  Eqs.  (12)  and 

(13). 

The  moving-source  data  along  the  track  shown  in  Fig.  1 
are  processed  using  the  above  method,  and  the  result  is 
shown  in  Fig.  11(b).  One  finds  zero  BER  except  near  the  end 
of  the  run  when  the  source-receiver  range  is  =£1  km.  The 
reason  is  attributed  to  the  fact  that  at  close  ranges,  the  vari¬ 
ous  multipaths  arrive  with  significantly  different  D/E  angles 
and  different  Doppler  shifts;  the  spread  of  the  Doppler  shifts 
creates  an  effective  Doppler  spread.  This  finding  was  sup¬ 
ported  by  Fig.  3(b),  which  indicates  a  decreased  channel  co¬ 
herence  time  or  a  higher  Doppler  spread  at  close  ranges.  To 
improve  the  BER  at  close  (=£  1  km)  ranges,  one  needs  to 
estimate  the  Doppler  shift  for  each  path,  which  is  signifi¬ 
cantly  more  difficult. 

D.  BER  as  a  function  of  decreasing  input  SNR 

To  assess  the  performance  of  the  transition  filter  at  a  low 
input  SNR,  ambient  noise  data  collected  during  the  same 
experiment  are  added  to  the  signal  data  (see  Ref.  10).  2  h  of 
noise  data  were  also  collected  producing  278  packets  of 
noise.  Multiple  packets  of  noise  data,  randomly  chosen,  were 
added  to  the  signal  data  to  generate  packets  of  data  with 
different  input  SNRs  ranging  from  -13  to  +30  dB.  The 
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FIG.  1 1 .  (Color  online)  Source-receiver  range  as  a  function  of  time  (packet 
number)  repeated  from  Fig.  3(a)  (upper  figure).  BER  as  a  function  of  packet 
number  (lower  figure). 


close  range  (=£  I  km)  data  are  excluded  from  this  analysis  as 
remarked  above.  BER  is  then  deduced  as  a  function  of  input 
SNR.  The  result  is  shown  in  Fig.  12(a)  as  a  function  of  the 
packet  input  SNR.  The  average  BER  is  shown  in  Fig.  12(b) 
as  a  function  of  input  SNR.  Compared  with  the  average  BER 
determined  from  the  fixed-source  data  [repeated  from  Fig. 
9(b),  shown  by  +],  one  finds  a  performance  degradation  of 
1-2  dB.  Note  that  while  the  cross-correlation  method  failed 
for  the  moving-source  data,  the  transition  filter  method 
works  pretty  well  for  both  fixed-  and  moving-source  data. 

IV.  PROBABILITY  OF  DETECTION  AND 
COUNTERDETECTION  RANGE 

The  probability  of  detection  ( PD )  depends  on  the  con¬ 
struction  of  the  detector;  the  more  the  detector  knows  about 
the  signal,  the  more  features  the  detector  can  use  to  improve 


the  detection.  If  the  signal  waveform  is  known,  one  can  em¬ 
ploy  a  coherent  detector,  e.g.,  a  matched  filter,  which  often 
yields  a  high  PD  given  the  processing  gain  derivable  from 
the  signal.  Noiselike  signals  are  difficult  to  detect  using  a 
coherent  detector  as  the  signal  is  random  and  difficult  to 
replicate  compared  with,  for  example,  an  impulselike  signal. 
This  also  favors  the  DSSS  signaling  as  compared  with,  for 
example,  frequency-hopping  frequency-shift-keying  signal¬ 
ing.  For  the  latter,  while  the  average  SNR  is  low,  the  SNR  for 
information  carrying  frequency  bins  is  high.  As  a  result,  the 
repeated  (or  time-varying)  frequency  hopping  pattern  may 
alert  the  operator  to  the  presence  of  a  man-made  signal.  Once 
the  operator  is  alerted,  more  sophisticated  signal  processing 
can  be  used  to  search  for  the  signal  parameters  and  construct 
an  improved  detector. 

If  the  signal  waveform  is  not  known,  one  is  left  with  an 
energy  detector.  A  commonly  used  energy  detector  is  the 
spectrogram,  which  detects  a  sudden  increase  in  the  SNR  as 
a  function  of  frequency  and  time.  This  normally  requires  a 
fairly  good  SNR  (say  5  dB)  over  the  signal  frequency  band 
and  is  most  effectively  done  by  a  human  being  who  can 
interactively  adjust  the  Fourier  transform  window  to  match 
the  signal.  A  computer-aided  detector  needs  to  know  the  sig¬ 
nal  bandwidth  and  signal  duration,  as  mismatch  in  bandwidth 
and  duration  between  the  signal  and  detector  can  signifi¬ 
cantly  decrease  the  detection  performance.  The  mismatch  in 
bandwidth  is  self-evident,  as  incorrect  bandwidth  results  in  a 
loss  of  part  of  the  signal  energy  and  addition  of  more  noise 
energy.  Likewise,  the  signal  duration  also  influences  the  de¬ 
tector  performance  (see  Sec.  IV  C).  If  the  integration  time  is 
much  longer  than  the  signal  duration,  it  effectively  reduces 
the  SNR  as  more  noise  is  included  in  the  detector  output. 
Without  knowing  the  signal  duration,  a  narrowband  detector 
is  often  used  to  detect  narrowband  energy  and  a  broadband 
pulse  detector  (i.e.,  short  integration  time)  is  commonly  used 
to  detector  a  burst  energy.  The  detectors  will  alert  the  opera- 


Input  SNR  (dB) 


Input  SNR  (dB) 


FIG.  12.  (Color  online)  BER  (a)  and  average  BER  (b)  as  a  function  of  input  SNR  using  the  transition  detector  method  applied  to  the  moving-source  data.  + 
denotes  the  fixed-source  data  of  Fig.  9(b).  The  solid  line  is  a  model  prediction  shown  in  Fig.  8(a).  Zero  BER  is  represented  by  10-6  on  the  logarithmic  scale. 
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tor  for  a  potential  signal  of  interest,  when  PD  exceeds  the 
recognition  differential  for  the  signal  of  interest. 

The  PD  for  an  energy  detector  depends  very  much  on  the 
input  SNR.  Pseudorandom  noiselike  DSSS  signals  with  an 
input  level  lower  than  the  ambient  noise  (e.g.,  SNR 
— 8  dB  within  the  signal  band)  are  difficult  to  detect  by  a 
narrowband  or  broadband  energy  detector  assuming  an  una¬ 
lerted  listener.  As  input  SNR  increases  (as  when  source- 
receiver  ranges  decrease),  the  probability  that  the  signal  will 
be  detected  increases.  From  a  system  point  of  view,  the  ques¬ 
tion  is  as  follows:  At  what  range  will  the  signal  be  detected 
or  what  is  the  counterdetection  range?  The  answer  is  channel 
dependent.  It  requires  a  determination  of  the  PD  and  the 
probability  of  false  alarm  (_PFA)  as  a  function  of  input  SNR, 
for  a  given  signal-fading  distribution.  This  topic  is  addressed 
in  Sec.  IV  A.  To  determine  the  input  SNR,  one  needs  to 
calculate  the  signal  transmission  loss  (TL)  as  a  function  of 
source-receiver  range  for  a  given  source  power  and  noise 
level.  This  calculation  is  done  for  a  typical  shallow  water 
environment  in  Sec.  IV  B.  This  section  also  shows  the  PD  as 
a  function  of  source-receiver  range  for  a  given  PFA.  The 
analysis  involves  two  receivers.  One  is  the  intended 
(friendly)  receiver  and  the  other  is  the  unfriendly  receiver. 


the  interceptor.  The  objective  is  to  minimize  PD  of  the  signal 
by  the  interceptor  for  a  given  PFA,  which  means  the  lower 
the  SL,  the  better  the  result.  On  the  other  hand,  the  SL  of  the 
transmitting  node  should  be  high  enough  so  that  the  message 
is  received  at  the  intended  receiver  with  a  high  enough  SNR 
for  decoding  the  message.  The  SL  is  a  system  parameter  that 
needs  to  be  determined  first.  See  Sec.  IV  B. 

A.  PD  and  PFA  of  an  energy  detector 

The  PD  and  PFA  for  a  broadband  pulse  energy  detector 
are  well  known  for  a  nonfading  channel15'16  and  a  Rayleigh- 
fading  channel,1  but  not  for  a  log-normal-fading  channel. 
The  data  as  discussed  above  indicate  that  the  fading  statistics 
for  the  DSSS  signal  in  an  underwater  acoustic  channel  fol¬ 
lows  the  log-normal  distribution.10  A  derivation  of  the  PD 
and  the  PFA  is  reviewed  below,  which  can  be  generalized  to 
include  the  log-normal  distribution. 

1.  Non  fading  channel 

For  a  nonfading  channel,  the  signal  (each  symbol)  has  a 
constant  amplitude  a.  Assuming  additive  white  Gaussian 
noise,  with  a  Rayleigh  amplitude  distribution,  the  probability 
of  detection  is  given  by  the  following  two  statistics: 


Pr(x)  : 


exp(-  (x2  +  fl2)/2a^)/0(ox/a^)  with  signal  present 


exp(-  x2/2a^,)  with  signal  absent. 


(14) 


where  crN  is  the  noise  variance  and  I0  is  the  modified  Bessel  function  of  zeroth  order.  Let  Y=X2  be  the  square-law  output  of 
an  energy  detector.  It  can  be  shown  that  the  probability  density  function  of  Y  becomes 


Priy) : 


— 7  exp(-  (y/2o^+  y))I0(2y]yy/2ofy  with  signal  present 

2  at 


(15) 


2  a) 


exp(-  yl2<j2N)  with  signal  absent. 


where  y=crl2a2N  is  the  mean  SNR  at  the  input  of  an  energy 
detector  (mean  input  SNR).  Given  a  threshold  rj,  let  the  PFA 
be  the  probability  that  the  noise  energy  is  higher  than  rj. 
Using  the  lower  equation  of  Eq.  (15),  one  has 
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PD=  —j  exp(-  {y/2a2N+  y))/0(2\/ yy/2o2N)dy 

J  v  2o^ 

r  oo 

=  exp(-  (u  +  y))I0(2\[yu)du,  (17) 

J  -ln(PFA) 

where  u=y  / 2o^  is  the  integration  variable  and  y  is  the  mean 
input  SNR. 


The  probability  of  detection,  PD,  is  the  probability  that  the 
signal  plus  noise  energy  is  higher  than  rj.  Using  the  upper 
equation  of  Eq.  (15),  it  can  be  expressed  as  a  function  of  the 
PFA  as  follows: 


2.  Fading  channel 

The  PFA  in  a  fading  channel  is  the  same  as  that  in  a 
nonfading  channel  [Eq.  (16)]  for  a  given  threshold  77.  The  PD 
can  be  expressed  as 
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r  oo 

PD  =  PD\sPs(s)ds 
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=  Ps(s)  I  — T  exp[-  (y  +  s2)/2a2N]IQ(s^y/a2N)dyds 

Jo  J  V 


—jexp{-y/2(T-N)A{y)dy, 
J  y  l(JN 


1 


(18) 


where  P D  s  is  the  conditional  probability  of  detection  for  a 
signal  with  an  instantaneous  signal  amplitude  s  as  given  in 
Eq.  (17),  ps  is  the  amplitude  distribution  of  the  signal,  and 

/•oo 

My)  =  exp(-  s1/2(j^)I0(syy/a^)ps(s)ds.  (19) 

Jo 


For  a  Rayleigh-fading  channel, 
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One  finds 
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where  77  is  the  threshold  set  by  PFA.  The  relationship  be¬ 
tween  PD  and  P¥A  can  be  expressed  as  follows: 

PFA=(PD)l+y,  (23) 

where  y=2cr1sl2<j2N=crlsl 0^  is  the  mean  input  SNR. 

For  a  log-normal-fading  channel,  the  amplitude-fading 
statistics  are  given  by  Eq.  (2)  for  the  signal  present  case. 
Substituting  the  log-normal  distribution  into  the  A(y),  Eq. 
(19),  one  finds 
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FIG.  13.  (Color  online)  Receiver  operation  characteristic  curves  showing  PD 
as  a  function  of  PFA  for  various  SNRs  and  signal-fading  distributions. 


The  relationship  between  the  detection  probability  and  the 
false  alarm  probability  is 


where  y=  E(S2) 1 2a^=  e2tls+2a's 1 2a^=  ^  2/2cr^is  the  mean  in¬ 
put  SNR  and 

Using  Eqs.  (17),  (23),  and  (25),  one  can  calculate  PD  as 
a  function  of  PFA  for  a  given  input  SNR,  the  so  called  re¬ 
ceiver  operation  characteristic  (ROC)  curve,  for  nonfading, 
Rayleigh-fading,  and  log-normal-fading  cases.  The  results 
are  shown  in  Fig.  13.  For  the  log-normal  fading,  we  have 
used  fj,s=-l. 6449  and  (7^=0.0304  obtained  from  the  distri¬ 
bution  of  input  signal  amplitude. 

The  ROC  curve  (Fig.  13)  is  commonly  used  for  detec¬ 
tion  analysis.  It  shows  the  dependence  of  PD  as  a  function  of 
PFA  for  a  given  SNR.  The  same  results  can  also  be  displayed 
in  a  different  manner,  for  example,  showing  the  PD  as  a 
function  of  SNR  for  a  given  PFA  as  will  be  done  in  the  next 
section  to  calculate  the  counterdetection.  One  notes  from 
Fig.  13  that  the  PD  versus  PFA  relation  based  on  the  log¬ 
normal  statistics  is  much  closer  to  the  nonfading  case  than  to 
the  Rayleigh-fading  case.  There  are  differences  between  the 
log-normal-  and  nonfading  cases,  which  are  not  clearly  seen 
in  a  log-log  plot. 


B.  Counterdetection  range 

As  PD  is  a  function  of  the  (mean)  input  SNR  and  the 
input  SNR  is  a  function  of  source-receiver  range  for  a  given 
SF,  one  can  calculate  the  PD  as  a  function  of  source-receiver 
range  for  a  given  PFA.  Counterdetection  range  will  be  de¬ 
fined  as  the  range  below  which  PD> 0.5.  To  model  the  SNR 
at  the  receiver,  we  calculate  the  TF  in  the  TREX04  environ¬ 
ment  using  the  range-dependent  acoustic  model  parabolic 
equation  (RAM  PE)  model.17  The  result  is  shown  in  Fig.  14 
for  a  17  kHz  signal.  Superimposed  is  a  sonar-equation-based 
calculation,  which  assumes  spherical  spreading  to  a  range  of 
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RANGE  (M) 


FIG.  14.  (Color  online)  TL  as  a  function  of  range  using  the  RAM  PE  model. 
The  smooth  line  is  the  empirical  model  calculation  based  on  sound  attenu¬ 
ation  by  sea  water. 

50  m  and  cylindrical  spreading  beyond  50  m.  For  absorption 
by  the  water  medium  we  use  the  following  expression  due  to 
Thorp:18 

/  0.1  f  40 f  ,,  \ 

a=[  — —  + - - — y  +  2.75  X  lCT4/2  + 0.003  X  0.875, 

\1  +f  4100  +f  ! 

(26) 

where  a  is  the  absorption  loss  in  dB/km,  /  is  the  acoustic 
frequency  in  kHz,  and  0.875  is  a  depth  correction  factor  de¬ 
termined  by  fitting  the  RAM  calculation;  the  original  for¬ 
mula  is  for  a  depth  of  1000  m.  To  determine  the  SNR  at  the 
receiver,  we  note  that  at  17  kHz,  the  noise  is  dominated  by 
wind  generated  noise,  with  a  spectral  level  given  by19 

10  log  N(f)  =  50  +  7.5w1/2  +  20  log/-  40  log (f+  0.4), 

(27) 

where  w  is  the  wind  speed  in  m/s.  Assuming  vv=  10  m/s,  one 
finds  NL=49  dB.  Given  the  TL  and  noise  level,  one  can 


Pfa  =0.01  Pd  Rayleigh 


Range  (km) 


calculate  the  SNR  as  a  function  of  range  for  a  given  SL. 
Given  the  SNR  one  can  then  calculate  the  PD  as  a  function  of 
range  for  a  given  PFA  using  the  equation  in  Sec.  IV  A.  The 
following  are  the  numerical  results  shown  in  Fig.  15  for  PD 
as  a  function  of  range  assuming  three  SLs:  143,  155,  and 
164  dB  and  Ffa=0.01.  These  SLs  are  needed  for  communi¬ 
cations  to  an  intended  (friendly)  receiver  at  a  range  of  ~2,  4, 
and  7  km.  Based  on  Fig.  12(b),  a  minimum  input  SNR  of 
-10  dB  is  required  at  the  intended  receivers  to  achieve  a 
minimum  (<1%)  BER.  The  minimum  input  SNR  is  used 
with  TL  to  determine  the  minimum  SLs. 

Figures  15(a)  and  15(b)  show  the  PD  as  a  function  of  the 
source-to-interceptor  range  for  Rayleigh  signal-fading  and 
log-normal  signal-fading  statistics,  respectively.  One  finds 
that  the  counterdetection  range  (for  PD^ 0.5)  is  approxi¬ 
mately  1.3,  3.6,  and  5.8  km  for  the  three  SLs  assuming 
Rayleigh-fading  statistics  and  approximately  1.4,  3.8,  and 
6.1  km  for  the  three  SLs  assuming  log-normal-fading  statis¬ 
tics.  Naturally,  the  higher  the  SL,  the  greater  the  counterde¬ 
tection  range.  As  expected,  the  counterdetection  range  for  an 
energy  detector  is  shorter  than  the  communication  range  to  a 
friendly  receiver  using  a  matched  filter.  Hence,  the  commu¬ 
nication  signal  will  not  likely  be  detected  by  an  interceptor 
located  outside  the  communication  range  (or  the  operation 
area).  But  as  the  interceptor  approaches  one  of  the  transmit¬ 
ting  nodes,  detection  by  the  interceptor  is  unavoidable.  The 
hope  is  that  the  interceptor  remains  unalerted,  such  as  when 
the  signal  is  noiselike. 

Note  that  at  PD= 0.5  the  differences  in  the  counterdetec¬ 
tion  ranges  are  relatively  small  between  the  two  signal¬ 
fading  scenarios.  However,  if  the  detection  ranges  were  de¬ 
fined  at  a  higher  PD,  e.g.,  PD=  0.9,  the  differences  between 
the  two  cases  would  be  significant.  The  difference  is  that  the 
PD  increases  faster  with  a  decreasing  range  for  the  log¬ 
normal  than  for  the  Rayleigh  distribution. 


Pfa  =0.01  Pd  lognormal 


Range  (km) 


FIG.  15.  (Color  online)  PD  as  a  function  of  source-receiver  range  for  PFA=0.01  for  three  assumed  SLs.  The  left  figure  assumes  Rayleigh  signal  fading  and 
the  right  figure  assumes  log-normal  signal  fading. 


J.  Acoust.  Soc.  Am.,  Vol.  124,  No.  6,  December  2008 


Yang  and  Yang:  Low  probability  of  detection  underwater  acoustic  3643 


C.  Integrating  over  multiple  observations 


If  the  signal  is  present  for  a  long  period  of  time,  the 
interceptor,  if  alerted,  can  improve  its  detection  by  integrat¬ 
ing  the  signal  energy  over  multiple  observations  (normally  in 
postprocessing  after  the  operator  has  been  alerted.)  Let  T  be 
the  sample  interval  and  y(t )  be  the  energy  at  time  t.  Assum¬ 
ing  that  all  samples  are  independent,  then  the  integrated  en¬ 
ergy  over  a  dwell  time  rd  can  be  approximated  as  follows: 

l  frd  ~ 1 

z=-|  y(t)dt~  2  y(kT),  (28) 

1  j  0  k=0 


where  Ns=rd/T  is  the  number  of  samples. 

For  large  Ns,  Z  is  approximately  Gaussian  distributed.  In 
contrast  with  Eq.  (16),  the  PFA  for  a  given  threshold  rj.  is 
given  by 
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where  /j,2=2Nsa 2N  and  o^=ANscrAN  are  the  mean  and  variance 
of  the  noise  only  case. 

Assuming  that  the  signal,  y,  follows  a  log-normal  distri¬ 
bution  (Sec.  IV  A),  then  the  probability  distribution  function 
of  z  is  given  by  (details  given  in  the  Appendix) 


P(z)  = 
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V27TCT, 
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(z  -  P-f 

2<L  . 


(30) 


where 

Pz  =  2Ns(jn(\  +  y), 

a{  =  4Nsa^l  +  2 y)  +  Ns  exp(4/rs  +  4cr^)[exp(4(7^)  -  1], 

(31) 

where  y=Es(s2) / 2cr^=exp(2fjLs+2a^) I2crh  is  the  mean  input 
SNR,  and  /jls  and  (r2  are  the  parameters  of  the  log-normal 
distribution  [Eq.  (2)].  In  contrast  with  Eq.  (25),  the  PD  is 
now  given  by  [see  Eq.  (All)] 

P  nl  Vz-2Nsofcl  +  y)  \ 

\  \4.Y,<f'|  1  +  2 y  +  y2(exp(4fr^)  -  1)]  / 

Using  Eqs.  (29)  and  (32),  one  plots  the  PD  as  a  function 
of  input  SNR  for  a  PFA=0.01  as  shown  in  Fig.  16  for  various 
values  of  Ns.  Comparing  Ns= 25  with  Ns=  1,  one  finds 
—  10  dB  improvement  in  SNR  at  PD= 0.5.  Thus,  knowing 
that  the  signal  bandwidth  and  hence  T  and  that  the  signal  is 
at  least  NST  long,  the  signal  could  be  detected  by  long  inte¬ 
gration  time.  As  an  example,  an  input  SNR  of  -8  dB,  which 
has  a  low  probability  being  detected  by  a  broadband  pulse 
detector,  could  now  be  detected  because  the  integration 
yields  an  effective  +2  dB  SNR.  This  assumes  that  the  opera¬ 
tor  has  been  alerted  that  there  is  a  man-made  signal  and  turns 
on  the  search  for  the  signal  parameters  and  long  integration 
detector. 


Input  SNR  [dB] 


FIG.  16.  (Color  online)  PD  as  a  function  of  input  SNR  for  PFA=0.01  for 
different  values  of  Ns. 

Figure  17  shows  the  probability  density  distributions  of 
the  noise  only  (upper  figure)  and  signal  plus  noise  (lower 
figure)  for  the  Ns=l  and  Ns=  25  cases.  For  the  Ns=  1  case, 
z=y,  one  finds  that  the  noise  density  distribution  decreases 
exponentially  with  energy.  With  the  signal  present,  while  the 
signal  (amplitude)  has  a  log-normal  distribution,  the  prob¬ 
ability  density  distribution  of  the  signal  plus  noise  decreases 
with  the  signal  +  noise  energy.  For  Ns=  25,  the  probability 
density  is  approximately  normally  distributed.  While  y  and 
z/Ns  have  the  same  mean  and  same  SNR  in  Fig.  17,  they 
have  different  probability  distributions.  The  fact  that  z  has 
the  same  probability  distribution  for  the  signal-plus-noise 
and  noise  only  cases  is  the  reason  that  PD  for  N .=25  is  large 
for  Ffa=0.01  (higher  than  that  for  Ns=  1,  for  the  same  PFA). 
One  may  argue  that  a  smaller  PFA  (e.g.,  PFA=0. 01/25  or  4 
X  10-4)  should  be  used  for  the  Ns= 25  case,  since  there  are 
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FIG.  17.  (Color  online)  Probability  distribution  of  noise  only  (upper  figure) 
and  signal  plus  noise  (lower  figure).  The  solid  line  is  for  Ns=\,  showing  the 
probability  as  a  function  of  y.  The  dashed  line  is  for  Ns=25,  showing  the 
probability  p(z)X^  as  a  function  of  z/Ns. 
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more  counts  of  false  alarms  for  Ns=25  than  for  Ns=l.  Note 
that  for  the  Ns= 25  case,  the  detector  does  not  require  the 
signal  to  be  detectable  for  each  sample  interval  T  but  only  for 
the  entire  time  period  of  25  T.  What  values  of  PFA  to  use  are 
system  dependent  and  likely  also  environment  dependent;  it 
is  beyond  the  scope  of  this  paper.  This  section  presents  only 
an  example.  One  can  calculate  PD  versus  PFA  based  on  the 
equations  in  the  Appendix. 

V.  SUMMARY  AND  DISCUSSION 

The  rapid  symbol  phase  change  as  a  function  of  time  for 
PSK  signals  has  been  well  noted  in  underwater  acoustic 
communications.  For  DSSS  signaling,  each  symbol  covers 
many  chips  and  the  symbol  phase  change,  accumulating  over 
the  chips,  can  be  large.  The  symbol  phase  change  can  be 
caused  by  a  random  medium,  which  produces  different 
phases  for  different  multipath  arrivals.  The  TREX04  data 
showed  that  the  differential  phase  between  adjacent  symbols 
for  fixed-source  and  fixed-receiver  transmissions  was  smaller 
than  ±90°,  suggesting  that  the  phase  imposed  by  the  random 
medium  changes  gradually  between  symbols.10  In  that  case, 
DPSK  signals  can  be  decoded  by  cross-correlating  the 
matched-filter  outputs  of  two  adjacent  symbols.  For  a  mov¬ 
ing  source  or  receiver,  the  fact  that  the  range  is  changing 
with  time  also  contributes  to  the  symbol  phase  change.  The 
additional  phase  change  contains  both  a  deterministic  as  well 
as  a  random  component;  the  latter  is  a  result  of  multipaths 
traveling  through  different  parts  of  the  random  medium  at 
different  times.  The  TREX04  data  showed  that  the  differen¬ 
tial  phase  between  adjacent  symbols  for  a  moving-source 
and  a  fixed-receiver  transmission  is  often  larger  than  ±90°. 
In  this  case,  the  cross-correlation  method  does  not  work  for 
the  moving-source  data. 

A  pair  of  transition  detectors,  based  on  the  same  prin¬ 
ciple  as  the  matched  filter,  are  proposed  in  this  paper  for 
moving,  as  well  as  for  fixed  source  data.  Whether  adjacent 
symbols  are  of  the  same  kind  or  opposite  kind  is  determined 
by  the  energy  of  the  transition  detectors.  The  transient 
(pulse)  detector,  being  an  energy  detector,  is  insensitive  to 
the  symbol  phase  fluctuations.  The  method  is  applied  to  the 
TREX04  data.  At-sea  noise  was  added  to  the  signal  data  to 
generate  data  of  low  input  SNR.  BER  was  measured  for  an 
input  SNR  from  +25  to  -15  dB.  The  moving-source  data 
yielded  a  comparable,  slightly  degraded  BER  compared  with 
the  fixed-source  data  using  the  same  method.  For  the  fixed- 
source  data,  the  performance  is  slightly  degraded  compared 
with  that  using  the  cross-correlation  method.  Whether  the 
performance  degradation  is  data  dependent  remains  to  be  in¬ 
vestigated. 

The  purpose  of  low  SNR  DSSS  communications  is  to 
minimize  the  probability  of  detection  and  interception  by  a 
hostile  interceptor.  Since  the  probability  of  detection  and  in¬ 
terception  is  detector  specific,  assuming  that  the  detector 
does  not  have  infinite  resources  to  employ  different  kinds  of 
detectors  simultaneously,  the  first  line  of  defense  is  to  mini¬ 
mize  the  probability  of  detection  by  the  two  most  popular 
detectors,  i.e.,  a  narrowband  energy  detector  and  a  broadband 
pulse  detector,  under  the  assumption  that  the  detector  has  no 


a  priori  knowledge  of  the  signal  parameters  or  the  signal 
parameters  are  constantly  changing  with  time.  Obviously, 
one  does  not  have  to  worry  about  narrowband  energy  detec¬ 
tor  for  DSSS  communications.  Anticipating  the  use  of  a 
broadband  pulse  detector,  the  question  is  at  what  range  will 
the  signal  be  detected,  the  counterdetection  range.  The  an¬ 
swer  is  case  dependent.  As  an  illustration,  we  assume  a  typi¬ 
cal  shallow  water  environment  and  analyze  the  probability  of 
detection  as  a  function  of  range.  Three  SLs  are  used  for  three 
intended  (friendly)  receivers.  Counterdetection  range  is  cal¬ 
culated  based  on  PD> 0.5  and  PFA=0.01.  The  analysis  as¬ 
sumes  a  single  receiver,  and  can  be  extended  to  an  array  of 
receivers  using  the  directivity  index.  As  expected,  the  lower 
the  SL,  the  shorter  the  counterdetection  range. 

If  the  interceptor  is  aware  of  a  man-made  signal  in  the 
water,  and  that  it  has  many  snapshots  of  data  to  investigate, 
the  interceptor  could  improve  its  detection  probability  by 
integrating  the  signal  energy  over  the  snapshots  of  data.  The 
probability  of  detection  using  a  long  integration  time  is  also 
analyzed  in  this  paper. 

Generally,  one  expects  the  lower  the  input  SNR,  the 
lower  the  probability  of  detection,  as  illustrated  by  the  nu¬ 
merical  examples  given  above.  But  communications  at  a 
lower  input  SNR  comes  with  the  price  of  a  lower  data  rate 
for  DSSS  signaling.  It  requires  a  longer  transmission  time  to 
get  the  same  throughput  compared  with  a  higher  data  rate, 
higher  SNR  communications  scheme.  One  expects  that  all 
signals  could  eventually  be  detected  if  the  detector  knows 
about  the  existence  of  the  signal  and  has  infinite  resources  to 
investigate  the  signal  properties.  By  way  of  illustration,  it  is 
pointed  out  that  low  SNR  DSSS  signal  could  be  detected  by 
long  integration  even  using  an  incoherent  detector,  if  such 
signals  are  known  to  be  present.  In  practice,  one  does  not 
have  infinite  resources.  The  advantage  of  low  SNR  DSSS 
communications  is,  from  the  practical  point  of  view,  that  it 
has  a  low  input  SNR  and  it  is  noiselike.  It  has  a  LPD  by  a 
broadband  pulse  energy  detector.  In  doing  so,  as  long  as  the 
operator  remains  unalerted  of  the  existence  of  the  signal,  the 
signaling  can  be  considered  covert. 
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APPENDIX:  PROBABILITY  OF  DETECTION  AND 
PROBABILITY  OF  FALSE  ALARM  WITH 
MULTIPLE  OBSERVATIONS 

Let  T  be  the  sample  interval  and  y(t)  be  the  energy  at 
time  t.  Assuming  that  all  samples  are  independent,  then  the 
integrated  energy  over  a  dwell  time  rd  can  be  approximated 
as  follows; 

\  rrd  N s~X 

z=-  y(t)dt~  2  y(kT),  (Al) 

1  j  0  k= 0 

where  Ns=rdIT  is  the  number  of  samples.  For  large  Ns,  Z 
may  be  approximately  Gaussian  distributed. 

Assuming  Gaussian  distribution  for  z,  one  has 
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where  Q  is  defined  in  Eq.  (29).  Thus  to  calculate  PD,  one 
needs  only  to  evaluate  the  mean  and  variance  of  the  signal 
plus  noise. 

With  the  signal  present,  the  mean  and  the  variance  of 
yk=y{kT)  for  a  given  instantaneous  signal  energy  s  in  the 
presence  of  AWGN  are  given  by 


/A 


■a  =  2  o^  +  s2  and  <ra  =  4o^(cr^+  s2). 


(A3) 


Letting  the  distribution  of  Yk  for  a  given  s  be  py  |j(>a|^'))  the 
distribution  of  Yk  becomes 


Pyfok)  =  J  Pyt\s(yk\S)Ps(S)ds- 
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Thus,  the  mean  and  the  variance  of  Yk  are  obtained  as  fol¬ 
lows: 
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For  a  Rayleigh-fading  channel,  using  Eq.  (20),  one  finds 
var(s2)=4cfj  and 


■Q 


V~2NscrN{\  +  y) 
^NsaUl  +  2y)  +  4  Nscr4 


(A9) 
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where  y=E/s2)l2crN=cr2sl  crN  is  the  mean  input  SNR.  The 
relationship  between  the  detection  probability  and  the  false 
alarm  probability  is  given  by 


Pd  =  Q 


Q-\Pfa)  -  \  J  Q-\Pfa)  -  yvW, 


\1  +  2y  +  Acr4IAcr4N 


=  Q 


1  +  y 
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For  a  log-normal-fading  channel,  using  Eq.  (2),  one 
finds  £'i(i2)  =  exp(2/u,s+2(7^)  and  var(s2)  =  exp(4/xs+4cr^) 
X[exp(4cr^)- 1],  Equation.  (A8)  then  yields 


Pd  =  Q 


y  -  2  A/J- 2n(  1  +  y) 


^Nsct4n[  1  +  2y+  y2(exp(4cr2)  -  1)]/ 


(All) 


where  y=£'J(s2)/2cr^=exp(2/u,J+2cr^)/2o^,  is  the  mean  input 
SNR.  The  relationship  between  the  detection  probability  and 
the  false  alarm  probability  is  given  by 

Q~1(Pfa)  -  yV V 


Pd  =  Q 
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and 
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where  E/s2)  and  var(,v2)  are  the  mean  and  variance  of  the 
signal  intensity,  and  y=E/s2)!2crN  is  the  mean  input  SNR. 
For  large  Ns,  assuming  a  Gaussian  distribution  for  z,  the 
mean  and  variance  of  z  are  as  follows: 


P\  =  2Nsa2N(\  +  y),  cr2  =  4Aicr^(l  +  2y)  +NS  var(s2) . 

(A7) 


The  corresponding  detection  probability  is 
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